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THE 100 MYR STAR FORMATION HISTORY OF NGC 5471 
FROM CLUSTER AND RESOLVED STELLAR PHOTOMETRY 
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We show that star formation in the giant Hn region NGC 5471 has been ongoing during the past 
100 Myr. Using HST/WFPC2 F547M and F675W, ground based JHKs, and GALEX FUV and NUV 
images, we have conducted a photometric study of the star formation history in the massive giant 
cxtragalactic Hn region NGC 5471 in M101. We perform a photometric study of the colour-magnitude 
diagram (CMD) of the resolved stars, and an integrated analysis of the main individual star forming 
clusters and of NGC 5471 as a whole. The integrated UV-optical-nIR photometry for the whole region 
provides two different reference ages, 8 Myr and 60 Myr, revealing a complex star formation history, 
clearly confirmed by the CMD resolved stellar photometry analysis. The spatial distribution of the 
stars shows that the star formation in NGC 5471 has proceeded along the whole region during, at 
least, the last 100 Myr. The current ionizing clusters are enclosed within a large bubble, which is 
likely to have been produced by the stars that formed in a major event ~20 Myr ago. 
Subject headings: H n Region, GEHR, star formation history, CMD 



1. INTRODUCTION 

Giant Extragalactic H n Regions (GEHRs) are impor- 
tant star formation tracers in the chemical evolution of 
galaxies and can be used as standard candles to estimate 
cosmological distances, thus they have been intensively 
studied. Equally important, they are the best indicators 
of the conditions that lead to massive star formation. 
These objects are very luminous Ha emitters. Their Ha 
luminosities, in the range 10 39 — 10 41 erg s~ 4 , imply large 
numbers of ionizing photons, 10 51 - 10 53 phs~\ an ion- 
izing power equivalent to up to several thousands of 05V 
stars. Such a massive stellar concentration provides an 
excellent laboratory to study the modes of massive star 
formation. 

The energy deposited in the interstellar medium in the 
form of stellar winds, supernova explosions and ultra- 
violet light produce complex stellar and gaseous spa- 
tial distributions which change rapidly during the few 
Myr following the onset of the first massive star gener- 
ation. Although the massive young clusters that power 
the GEHRs are expected to form on a timescale shorter 
than 10 Myr. some studies in regions such as 30 Doradus 
in the LMC reveal that their stars may not be coeval 
(Walborn ct al. 1999). The first population formed in 
the central regions may have induced star formation in 
the periphery through the action of UV radiation and 
stellar winds that compress the surrounding molecular 
clouds. Therefore, the study of the interrelation between 
stars, gas and dust components is highly relevant to un- 
derstand the nature of these regions. 

30 Doradus in the Large Magellanic Cloud (LMC) and 
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NGC 604 in M33 are the two largest GEHRs in the Lo- 
cal Group and they have been the subject of very de- 
tailed analyses. Both are dominated by very young (2-4 
Myr) ionizing stellar populations, but they present two 
rather different stellar distributions. While NGC 604 is 
an extended Scaled OB Association (SOBA) with ~200 
O+WR massive stars surrounded by nebular filamentary 
structures, 30 Doradus is powered by a Super Star Clus- 
ter (SSC) - a much more compact core (Mafz Apellaniz 
2001) - and a halo, with different stellar populations 
ranging from 2 to 20 Myr (Walborn & Blades 1997). 
Both GEHRs include an ongoing star formation with 
young stars still embedded in their parent molecular 
cloud. The structural similarities and differences be- 
tween these well-studied and spatially well-resolved Local 
Group GEHRs have made them obvious ideal laborato- 
ries in which to study massive star formation, but in 
order to extrapolate our knowledge to distant luminous 
unresolved starbursts it is critically necessary to go a step 
further and to analyze in detail the somewhat more dis- 
tant GEHRs just outside the neighbourhood of the Local 
Group. 



M101 (NGC 5457) is a giant spiral galaxy located at a 
distance of 7.2 Mpc (Stetson et al. 1998) yielding a lin- 
ear scale of 34.9 pc/arcsec. This galaxy contains a large 
number of very luminous GEHRs, of which NGC 5471 
(a 200 o = 14 7l 04 m 29 s , S 2Q0 q = +54°23'48"; 1 = 101?78, b 
= +59?63) is one of the outermost, at a galactocentric 
radius of about 25 kpc (see figure 1). The Ha morphol- 
ogy of NGC 5471 shows multiple cores with surrounding 
nebular filaments, extending over a diameter of ~ 17", or 
~ 600 pc. Ground-based photometry shows five bright 
knots, designated by Skillman (1985) as components A, 
B, C, D and E. Component A is as luminous as 30 Do- 
radus, while components B, C and E are comparable to 
NGC 604, the two prototypical spatially resolved mini- 
starbursts in the Local Group. Thus NGC 5471 is truly 
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Fig. 1. — Large scale image of M101. The grayscale is the Digital 
Sky Survey image and the contour image is the H I from the VLA 
(Braun 1995). NGC 5471 is located in the outskirts of the galaxy, 
in an arm concentration of H I. 

a giant star forming region. 

All these features make of NGC 5471 an excellent can- 
didate for the study of complex star formation, being 
close to the limit distance between the spatially well re- 
solved Local Group GEHRs and the more distant unre- 
solved starbursts, thus providing an intermediate step in 
our understanding of very massive star forming regions. 

This paper is organized as follows. Section 2 describes 
the observations and data reduction. Section 3 presents 
the integrated analysis of both the individual large clus- 
ters and the whole of NGC 5471, and the photometry 
of the resolved stars. The results found in section 3 arc 
elaborated and discussed in Section 4. Finally, the con- 
clusions are summarized in Section 5. 

2. OBSERVATIONS AND DATA REDUCTION 

We have retrieved ultraviolet NUV and FUV images 
of M101 from the GALEX archive, as well as broad and 
narrow band optical images for NGC 5471 from the HST 
archive. Finally, we have obtained near infrared JHK 
images. 
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Fig. 2.— HST WFPC2 images of NGC 5471 in continuum cor- 
rected nebular filters Ha (F656N), [SB] (F673N), and emission line 
corrected continuum filters F547M and F675W'. The five knots 
from Skillman (1985) are marked in the Ha image. An ellipse in- 
dicating the core of the region is drawn in the Ha and F675W' 
images. 

2.2. HST/WFPC2 Imaging 

Hubble Space Telescope (HST) Wide Field Planetary 
Camera 2 (WFPC2) images of NGC 5471 were retrieved 
from the HST archive. The images were taken on 1997 
November 1 through two emission-line filters: F656N 
(Ha) and F673N QS II]), and two continuum filters: 
F547M (Stromgren y) and F675W (WFPC2 R) (Heyer 
et al. 2004). Each camera images onto a Loral 800x800 
CCD which gives a plate scale of 0'.'046 pixel -1 for the 
PC camera and O'.'IO pixel -1 for the three WF cameras, 
with a readout noise of ~5 e" and a gain of 7 e~ /DN 
for this observations. NGC 5471 is centered in the WF3 
camera in all images. The scale of the WF CCDs at 
NGC 5471, for which we assume a distance modulus of 
(m - M) = 29.3 (Stetson et al. 1998), is 3.5 pc pixel" 1 . 
Table 1 lists the details concerning the WFPC2 data. 



2.1. GALEX data 

The Galaxy Evolution Explorer (GALEX) far- 
ultraviolet (FUV; X ref = 1530 A, AA = 400 A) and 
near-ultraviolet (NUV; \ ref = 2310 A, AA = 1000 A) 
images of M101 were retrieved from the Nearby Galaxies 
Survey (NGS). 

The GALEX instrument is described by Martin et al. 
(2005) and its on-orbit performance by Morrissey et al. 
(2005). GALEX FUV and NUV imaging was obtained 
with total exposure times of 1041 s in each band. The 
lcr NUV (FUV) sensitivity limit of the GALEX images 
is 27.6 (27.5) AB mag arcsec -2 for these data (Bianchi 
et al. 2005). At the distance of M101, the linear scale of 
34.9 pc/" implies that the GALEX point spread function 
(PSF) of 4.6 arcsec translates to 160 pc; this is barely suf- 
ficient to resolve some of the structure in NGC 5471, and 
we will only use the GALEX data to study the integrated 
photometry. 



The HST-pipelinc WFPC2 images were subjected to 
the usual processing using the IRAF 5 and STSDAS 6 soft- 
ware packages, following the procedures outlined in the 
HST Data Handbook. We used the multiple exposures 
to correct for the incidence of cosmic rays as well as to 
eliminate hot pixels and other defects on the CCD. We 
finally co-added the images to produce deep mosaics. 

Given the luminosity of NGC 5471, its Ha and [SII] 
emission contributes significantly to the broad band 
(F675W) image. To eliminate the nebular contribution 
from the broad band and to produce continuum-free Ha 

5 IRAF is distributed by the National Optical Observatory, 
which is operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 

6 Space Telescope Science Data Analysis System (STSDAS) is 
a product of the Space Telescope Science Institute, operated by 
AURA for NASA 



The continuous SFH of NGC 5471 



3 



TABLE 1 

Journal of HST/WFPC2 observations of NGC 5471. 



Filter Band Exposure (s) Observation ID 

F547M Stromgren y 2 X 600 U4DN030DR, U4DN030ER 

2 x 100 U4DN030FR, U4DN030GR 

20 U4DN030HM 

F675W WFPC2 R 2 x 400 U4DN0308R, U4DN0309R 

2 x 50 U4DN030AR, U4DN030BR 

10 U4DN030CR 

F656N Ha 180 U4DN0305R 

2 x 600 U4DN0306R, U4DN0307R 

F673N [S II] 3 x 700 U4DN0301R, U4DN0302R, U4DN0303R 

300 U4DN0304R 



Note. — All images were obtained on 1997 November 1 for the cycle 6 program 
GO-6928, with You-Hua Chu as PI. The object of study is centered in the WF3 



and [SII] emission line images, we followed the procedure 
decribed below. First we multiplied the images by the 
PHOTFLAM parameter obtained by synphot to convert 
from count rates to flux densities, and then divided them 
by their respective time exposure. From the SYNPHOT 
task in STSDAS we determined the filter widths, and 
multiplied the flux densities found by the corresponding 
filter widths in order to obtain fluxes. The correct filter 
widths are the filter rectangular width RECTW (Lurid- 
iana et al. 2003) defined as the width in Angstroms of a 
rectangle with the same total area as the total transmis- 
sion curve and a height given by the peak in the curve. 
We then scaled the Ha, [S II] and F675W images and per- 
formed an iterative subtraction to obtain a new, emission 
line free, broad band image (F675W). 

Figure 2 shows the Ha image. The five brightest com- 
ponents first noted by Skillman (1985) and designated as 
A, B, C, D and E, are distinctly observed in that figure. 
The presence of a supernova remnant (NGC 5471B) is 
clearly seen as an intense [SII] emission source in the B- 
component in the F673N image. Previous works have as- 
sociated X-ray sources with supernova remnants (SNRs) 
in M101. Chen et al. (2002) conclude that the energetic 
SNR in NGC 5471B (first identified by Skillman 1985) 
was very likely produced by a hypernova, an explosion 
one or two orders of magnitude brighter than the canon- 
ical 10 51 erg energy explosion associated with a normal 
SN, although this result has been questioned recently 
(Jenkins et al. 2004). 

2.3. Near-infrared photometry 

JHK S photometry was performed on 2000 April 29, 
using the 3.58m TNG (Telescopio Nazionale Galileo) 
at the Observatorio del Roque de los Muchachos in 
La Palma. The camera used was ARNICA (ARcetri 
Near Infrared CAmera), with a CCD detector NICMOS 
256x256 HgCdTe and gain of 19 e"/ADU. The spa- 
tial scale was 0'.'35 pixel -1 and the field of view 90x90 
arcsec 2 . The details of these observations are given in 
Table 2. 

The flux calibration of the images has been done using 
2MASS images of NGC 5471, retrieved from the 2MASS 
Extended Source Catalog. We first measured the inte- 
grated flux of the 2MASS images in a circle contain- 
ing the region and then, we used that aperture centered 
in the same coordinates and covering the same area in 



TABLE 2 

Journal of TNG/ARNICA observations of 
NGC 5471 



Band 


Time (UT) 


Duration (s) 


air mass 


J 


04:09:00 


1 X 60 


1.38 


11 


01:03:08 


3 x 15 


1.11 


K s 


21:59:29 


10 x 5 


1.31 



the ARNICA images. We obtained the photometric con- 
stants comparing the calibrated 2MASS fluxes with those 
from ARNICA by means of the transformation equations 
of the two systems. 

In order to compare the two sets of optical and NIR 
data, we aligned the HST and TNG images. We de- 
graded the WFPC2 images to the spatial scale and seeing 
of the IR ones, then we identified several objects in both 
data sets and marked them with the IRAF center rou- 
tine. Next, we aligned the images with the geomap and 
geotran routines, using the HST images as reference. 



Figure 3 shows the images and contour plots of the in- 
frared bands. The multiple knots present in these images 
are well identified with those in Figure 2 in the F547M 
and F673W filters. 

3. RESULTS 

3.1. Resolved Stellar Photometry 

3.1.1. WFPC2 Photometry 

The stellar photometric analysis was performed with 
the HSTphot package (Dolphin 2000a). This package is 
specifically designed for use with HST WFPC2 images. 
It uses a library of Tiny Tim (Krist 1995) undersampled 
point-spread functions (PSFs) for different locations of 
the star on the camera and of the star within the pixel, 
to center the star and to find its magnitude, given in the 
flight system magnitude. It also contains several pro- 
grams for general reduction procedures. The first step 
is to run the mask task to mask out the bad pixels and 
other image defects identified in the data quality images. 
The sky computation is made by getsky, which takes all 
pixels in an annulus around each pixel, determines the 
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Fig. 3. — JHK S grayscale images and contour plots. The isophotes 
ellipse indicating the core of the region is drawn in the H image. 

sky value and calculates the sky background map. The 
next step is to run crmask for cosmic ray removal, which 
uses the sky file created by getsky. It has the capabil- 
ity of cleaning images that are not perfectly aligned, and 
it can handle images from different filters. To combine 
multiple images per filter and to produce a final deep 
image per filter, the coadd procedure was used. The fi- 
nal step requires the hotpixels procedure, which uses 
the result from getsky and tries to locate and remove 
all hot pixels. This is an important step, since hot pixels 
can create false detections and also, can throw off the 
PSF solutions. 

The fistpfiot routine was run on the images in the 
F547M and F675W' bands. This task performs stel- 
lar PSF photometry on multiple images from different 
filters, including alignment and aperture corrections, as 
well as PSF modifications to correct for errors of geomet- 
ric distortion via the Holtzman et al. (1995) distorsion 
correction equations and the 34th row error, noted by 
Shaklan et al. (1995) (see also Anderson & King 1999), 
and correction for charge transfer inefficiency (Dolphin 
2000). There are several parameters to use during the 
run of hstpriot. We enabled the determination of the 
"local sky" (from now on, nebular case) , because of the 
rapidly- varying background of NGC 5471, but we also 
used, in a different run, the usual (stellar) fit done by 
getsky. We made a cross-correlation between the nebu- 
lar and stellar lists of stars in order to match them and 
find the sources detected using both types of sky deter- 
mination. To have a better estimate of the sky for those 



range from 20 up to 100 per cent of the peak value of the image. An 

matched stars, we took as magnitude the mean of both 
measurements m = (m ne i, u i ar + m s teiiar)/2 and as error: 

/ \m n -ms\\ 
a m = max I a mn , a ms , I 

The final list of stars is made up of: i) matched stars; ii) 
objects from the stellar list not included in the matched 
star list (best detected in less nebulous regions); and iii) 
stars from the nebular list with no companion in the 
stellar list. 

We selected "good stars" from the hstphot output. 
Object types were classified as good star, possible un- 
resolved binary, bad star, single-pixel cosmic ray or hot 
pixel, and extended object. Less than 10% of the stars 
were not classified as "good stars". To ensure that we 
only have point-like sources, we also used the "sharp- 
ness" parameter and a measure of the quality of the fit 
(X 2 ) to reject false star detections in regions with struc- 
tured nebulosity or artifacts. The final number of stars 
detected in the WF3 chip with these parameters in both 
filters were 1261. 

3.1.2. Selection Effects 

Completeness — In order to face the study of the star for- 
mation history (SFH) of an object and simulate synthetic 
diagrams, it is compulsory to undertake a completeness 
analysis. A thorough description of the completeness 
analysis can be found in Stetson & Harris (1988) and 
Aparicio & Gallart (1995). We briefly describe below its 
basis and main steps. 
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Incompleteness is mainly due to two factors: overlap- 
ping of stellar profiles in crowded regions and background 
noise. These effects are more important for faint stars in 
dense stellar regions. Consequently, in a crowded cluster 
incompleteness will depend on apparent magnitude and 
spatial position. The best way to account for crowding 
effects is to add artificial stars to the images and to sub- 
ject them to the same sample selection as for real stars. 
The number of injected stars must be sufficiently high 
as to provide acceptable statistical accuracy, but not so 
high as to change the level of crowding. To avoid over- 
crowding effects, no more than some 10% of the number 
of real stars should be added as artificial stars. These 
new artificial images (synthetic frames) are processed in 
the same way as the real ones, and the magnitudes of 
all the stars are obtained following the same criteria and 
using the same parameters. At this point, we will have 
three different tables: i) the first one containing magni- 
tudes, positions and colors of the real stars (real stars); 
ii) the second one including all the data for the injected 
stars, which are known beforehand; Hi) the third list 
containing the data from the synthetic frames. 

The synthetic table is then cross correlated with the 
other two (real + injected), producing a new table (com- 
pleteness) containing the stars matched from the arti- 
ficial list. Therefore, in this final list we have the ini- 
tial and recovered magnitudes and colors of the injected 
stars, as well as their coordinates. The table will also 
contain stars which have not been paired, and a flag in- 
dicating that they have been lost. The ratio between the 
number of simulated stars recovered and injected pro- 
vides the estimate of the completeness factor for each 
magnitude bin. 

A grid of artificial stars was generated on a 2- 
dimensional CMD and distributed according to the flux 
of the images with an artificial star routine provided by 
HSTphot. The parameters of the routine are the mini- 
mum and maximum of the measured color and magni- 
tude. The magnitude steps used were multiple of 0.5, 
while color steps are by 0.25. We asked the routine to 
distribute artificial stars on the CMD in accordance with 
the number observed. Approximately 40,000 fake stars 
were added in each image (at different trials, in order to 
leave the crowding conditions unaltered) and were given 
random magnitudes and colors in the observed range. 
We estimated the completeness in two different regimes 
of crowding: inside the core of the region (inner ellipse in 
Figure 4) and in the rest of the image. For the core, the 
50% completeness of the F547M filter luminosity func- 
tion derived on the basis of the CMD data is reached at 
24.8 magnitude, while for the F675W' filter is 24.3 mag. 
For the rest of the image, the 50% completeness of the 
F547M filter is reached at 25.2 magnitude, while for the 
F675W filter is 24.7 mag. 



Field star subtraction — Statistical decontamination has 
been widely used to accomplish the field star subtracion 
of CMDs. The idea is very simple: our region of study, 
NGC 5471, consists of the sum of the intrinsic popu- 
lation, IP, and the contaminating population, FS, and 
therefore our CMD shows both populations. The field 
star contribution can be estimated using the CMD of an 
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Fig. 4. — Continuum free Ha image of NGC 5471. The circular 
aperture of 15" where the integrated photometry was performed to 
obtain the total flux is overplottcd together with the inner ellipse 
defining the core of the region. 

adjacent field FS' containing only foreground and back- 
ground stars, assuming that it is statistically represen- 
tative of FS. Accordingly, the FS contribution may be 
removed by comparing the local density of stars in the 
two diagrams (Bellazzini et al. 1999): 

CMD(IP + FS) - CMD(FS') ~ CMD(IP) 

As NGC 5471 is centered in the WF3, we assume that 
all of the stars in the other chips are field stars. This 
is not completely true because this GEHR is in the out- 
skirts of an M101 arm, where a few traces of star forma- 
tion activity can still be found. Nevertheless, they are 
easily detected and no more than 20 stars in the WF2 
(the only chip used to determine the FS population) were 
found over the chip. The rest of the stars showed a ho- 
mogeneous distribution and were considered a sample of 
field stars. A CMD(IP+FS) was constructed for a re- 
gion inside a circular apcture of 15" which includes the 
core+halo structure of the region. 

To perfom the statistical decontamination, we followed 
the procedure described by Mighell et al. (1996) and Bel- 
lazzini et al. (1999). For a given star in the CMD(IP+FS) 
with coordinates [to ± o~ m ; c ± cr c ] , we count how many 
stars (Nip+fs) can be found within the ellipse of axes 
[max(a m , 0.1); max(a c , 0.1)] in the original diagram. We 
also counted the number of stars (Nps') found in the field 
CMD(FS') within the same ellipse centered at the same 
coordinates. 

The probability P that a given star in the CMD(IP + 
FS) is actually a member of NGC 5471 is defined as 
follows: 



P 



1 — vain ( — , 1.0 



N 



IP+FS 



where a is the ratio of the area of the IP+FS field (a 
circular aperture of 15") to the area of the FS' field (in 
this case, the whole WF2 chip). We can determine the 
probable membership of this star by comparing P with 
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Fig. 5. — Decontamination of the NGC 5471 observed CMD 
[panel (a)] using the field stars CMD [panel (b)]. The decontami- 
nated NGC 5471 CMD is shown in panel (c) (610 stars) while the 
CMD of the rejected stars is displayed in the panel (d) (73 stars). 

an uniform random number < P' < 1; if P' < P then 
the star is accepted as a member of NGC 5471, other- 
wise it is rejected and considered as a field star. Since 
this cleaning method is probabilistic in nature, a decon- 
taminating trial represents only one out of an infinite 
number of different realizations of the cleaned CMD for 
NGC 5471. 



Figure 5 shows the result of this process of decontam- 
ination for the CMD of NGC 5471. Panel (a) displays 
the CMD observed (with both intrinsic and field stars), 
panel (b) the field stars, panel (c) the decontaminated 
CMD of and panel (d) the CMD of the stars rejected. 
Table 3 shows the magnitudes of the first 20 stars of the 
resulting clean CMD of panel (c); the full table can be 
accessed in the electronic edition. 

3.1.3. Isochrones 

To study the CMD we use the isochrone set that cor- 
responds to the grid by Girardi et al. (2002) 7 , covering 
the age range 6.60 < logt < 10.25 and mass range 0.15 
to 100 M Q , based on the Z=0.004 (1/5 Z Q ) models (see 
end of section 3.2 for a discussion on the choice of the 
metallicity) from Girardi et al. (2000) and Bcrtclli et al. 
(1994) that include overshooting and a simple synthetic 
evolution of TP-AGB (Girardi and Bcrtclli 1998). Atmo- 
sphere models arc from BaSeL 3.1 WLBC 99 corrected 
library (Westera et al. 2002) at the same metallicity as 
the isochrone set. 

The procedure is as follows: (1) we obtained the mag- 
nitudes in the HST filters for all stellar models in the 
atmosphere grid. (2) For each point in the isochrone we 
obtained the stellar flux in each of the bands by a linear 

7 Available at http://pleiadi.pd.astro.it/ 
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Fig. 6. — Final CMD corrected for field star contamination, 
together with five isochrones from 4 to 70 Myr. The line across the 
bottom shows the completeness limit at the 50% level as derived 
from false star tests for the core (dashed line) and for the rest of 
the image (solid line). The right-hand side scale gives the absolute 
magnitudes corresponding to a distance of 7.2 Mpc. 

interpolation in the log g — logTcff grid. (3) The final 
luminosity was obtained by multiplying by the square 
radius of the star. For log g — logTeff values not covered 
by the atmosphere grid, the boundary values have been 
used with a caution flag. In terms of the present work, 
these stars do not modify the results obtained. 



In Figure 6 we show the final CMD (corrected for field 
star contamination) . The 50% completeness limit is also 
plotted, and five isochones from 4 to 70 Myr are shown 
in colour. The right hand side scale gives the absolute 
magnitude corresponding to a distance of 7.2 Mpc. 

3.2. Cluster Photometry 

We have derived the integrated magnitudes and col- 
ors in eleven different apertures covering the main com- 
ponents in the core of NGC 5471. The apertures were 
chosen and defined in the H-band image, maximizing the 
area in order to include the features of all images, since 
there were some small offsets in the position and inten- 
sity of the knots in the optical and infrared images re- 
spectively. Then, the same polygons were used for the 
J, K s and HST images. Polygons were marked by using 
the IRAF task polymark and the integrated magnitudes 
were computed using the IRAF task polyphot. The in- 
tegrated magnitudes of the HST images are given in the 
VEGAMAG system. The polygons defined are shown in 
Figure 7. 



The integrated magnitudes and colors derived in each 
area were corrected for the Galactic foreground extinc- 
tion toward M101 using the color excess E{B — V) = 0.01 
(Schlegel et al. 1998), although its effect is negligible. 
The absorption coefficients for the 5 bands were derived 
assuming Ry = Ay / E(B — V) = 3.1, and using the 
Aa curve from Cardelli et al. (1989). The uncertainty 
of the individual fluxes have been computed including 
the Poisson photon count error, the uncertainty on the 
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TABLE 3 
STELLAR PHOTOMETRY 



X 


Y 


R.A. 


Decl. 


m _F547M 


m F675W 


m F547M — m F675W 


ID pixel 


pixel 


(J2000) 


(J2000) 


(mag) 


(mag) 


(mag) 



1 


359 


56 


562 


57 


14 


04 


27 


50 


54 


23 


11 


53 


24 


21 


± 





07 


23 


45 


± 





05 


0.79 


± 


0.09 


2 


324 


75 


512 


37 


14 


04 


27 


56 


54 


23 


17 


59 


24 


57 


± 





10 


23 


93 


± 





07 


0.65 


± 


0.12 


3 


352 


11 


544 


56 


14 


04 


27 


57 


54 


23 


P: 


36 


24 


05 


± 





07 


23 


95 


± 





06 


0.09 


± 


0.09 


4 


332 


62 


514 


12 


14 


04 


27 


61 


54 


23 


46 


96 


21 


59 


± 





09 


21 


88 


± 





19 


-0.30 


± 


0.21 


5 


378 


56 


569 


56 


14 


04 


27 


62 


54 


23 


39 


79 


24 


04 


± 





06 


23 


62 


± 





06 


0.42 


± 


0.08 


6 


323 


33 


496 


14 


14 


04 


27 


66 


54 


23 


48 


91 


24 


07 


± 





06 


21 


13 


± 





12 


-0.06 


± 


0.14 


7 


373 


56 


553 


15 


14 


04 


27 


69 


54 


23 


41 


37 


21 


18 


± 





07 


23 


21 


± 





05 


0.94 


± 


0.09 


8 


293 


29 


455 


69 


14 


04 


27 


69 


54 


23 


53 


94 


23 


62 


± 





05 


22 


57 


± 





03 


1.05 


± 


0.06 


9 


351 


67 


523 


45 


14 


04 


27 


71 


54 


23 


45 


04 


23 


85 


± 





07 


23 


87 


± 





11 


-0.02 


± 


0.13 


10 


290 


75 


449 


59 


14 


04 


27 


71 


54 


23 


54 


57 


23 


86 


± 





06 


23 


02 


± 





01 


0.84 


± 


0.07 


11 


368 


80 


540 


15 


14 


04 


27 


74 


54 


23 


42 


65 


23 


14 


± 





03 


23 


21 


± 





01 


-0.10 


± 


0.05 


12 


346 


68 


513 


19 


14 


04 


27 


74 


54 


23 


16 


12 


25 


26 


± 





15 


21 


89 


± 





17 


0.38 


± 


0.23 


13 


381 


63 


555 


86 


14 


04 


27 


75 


54 


23 


40 


66 


21 


08 


± 





07 


21 


18 


± 





11 


-0.10 


± 


0.13 


11 


360 


90 


530 


61 


14 


04 


27 


75 


54 


23 


13 


91 


23 


96 


± 





06 


21 


09 


± 





09 


-0.14 


± 


0.10 


15 


318 


82 


469 


76 


14 


04 


27 


81 


54 


23 


51 


25 


25 


22 


± 





15 


25 


01 


± 





20 


0.18 


± 


0.26 


16 


312 


17 


461 


48 


14 


04 


27 


81 


54 


23 


52 


29 


23 


20 


± 





04 


23 


22 


± 





04 


-0.02 


± 


0.05 


17 


291 


51 


435 


95 


14 


04 


27 


81 


54 


23 


55 


57 


25 


38 


± 





21 


25 


23 


± 





19 


0.16 


± 


0.28 


18 


359 


87 


517 


56 


14 


04 


27 


83 


54 


23 


11 


98 


24 


20 


± 





09 


21 


26 


± 





10 


-0.06 


± 


0.13 


19 


354 


84 


509 


85 


14 


04 


27 


84 


54 


23 


15 


89 


21 


02 


± 





06 


23 


67 


± 





05 


0.34 


± 


0.08 


20 


319 


56 


460 


53 


14 


04 


27 


88 


54 


23 


51 


91 


24 


15 


± 





08 


21 


87 


± 





20 


-0.42 


± 


0.22 



Note. — Stellar photometry. Table 3 is published in its entirety in the electronic edition, 
portion is shown here for guidance regarding its form and content. 
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Fig. 7. — Plot of the eleven polygons, limiting the different ar- 
eas where the integrated magnitudes are computed. The numbers 
belong to our own definition while the A, B, C, D and E letters are 
the knots identified by Skillman (1985). The hypernova is located 
in knot 6. 

background setting and the photometric calibration un- 
certainty. The results of the integrated photometry are 
shown in Table 4. 

We have also measured the Ha flux for each poly- 
gon. To measure the integrated Ha flux of the whole 
region, we have decontaminated the F656N image from 
the broad band F675W contribution. We have then plot- 
ted the contour at 3er from the sky (with isophotal flux 
6.3 x 10 -15 erg cm" 2 s" 1 arcscc -1 ) and found the circular 
aperture centered in the region that best fits that isopho- 
tal contour. The radius of this aperture is 15 arcsec from 



the center (see Figure 4). The flux measured inside this 
circular aperture is (3.65 ± 0.17) x 10 -12 erg cm~ 2 s -1 . 
We have used the PHOTFLAM value given by synphot 
and a value of RECTW of 28.3 A. Measured Ha fluxes 
for each of the polygonal selections and the parameters 
derived from them are given in Table 5. The effective 
radius (R e /y) of a polygon is defined as the radius of a 
circular aperture of the same area. 

Since GALEX images do not have the necessary resolu- 
tion, we were not able to get the fluxes for each polygon. 
As in the case of the nebular flux, we measured the FUV 
and NUV fluxes inside 15" radii apertures centered in 
NGC 5471. We obtain a FUV flux of (9.9 ±0.1) x 10~ 14 
erg cm~ 2 s _1 A -1 (14.19 ± 0.05 AB magnitudes) and 
a NUV flux of (4.5 ± 0.1) x 10" 14 erg cm" 2 s" 1 A" 1 
(14.24 ± 0.05 AB magnitudes). 

In order to use these values in the analysis, we need 
to transform AB magnitudes to the VEGAMAG system. 
The transformation formula is (Mafz Apellaniz 2007): 



AB„ 



ZP p )+AB p (Vega) 



(1) 



where AB p is the magnitude of the object in the AB 
system and filter passband p, m p the magnitude of the 
object in filter p using Vega as reference spectral energy 
distribution (SED), ZP P is the zero point of the p filter 
in the Vega system and AB p (Vega) is the AB magnitude 
of Vega in filter p. This magnitude is given by: 



AB p (Vega) = -2.hlog 1 



(X)XdX\ 
J P p (X)(cf^ AB /X)dX J 



(2) 

where P P (A) is the dimensionless filter sensitivity func- 
tion of filter p, t\y e ga(X) is the Vega spectrum and I v .ab 
= 3.63079 x 10 -20 erg s^ 1 cm~ 2 Hz" 1 is the reference 
spectral energy distribution of the AB system. 
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TABLE 4 
INTEGRATED PHOTOMETRY 



Aperture 


J 


H 


K s 


m 547M 


m 675W 


m 675W 



N5471- 


1 


16.43 


± 


0. 


11 


16.26 


± 


0. 


12 


16.03 


± 





14 


16 


66 


± 





06 


15, 


56 


± 





06 


16. 


30 


± 





,06 


N5471- 


■2 


16.90 


± 


0. 


11 


16.57 


± 


0. 


12 


16.48 


± 





14 


17 


21 


± 





06 


16 


30 


± 





06 


16. 


94 


± 





.06 


N5471- 


■3 


17.35 


± 





11 


17.09 


± 


0. 


12 


17.12 


± 


0. 


14 


17 


22 


± 





06 


16 


81 


± 





06 


17 


11 


± 





,06 


N5471- 


■1 


17.52 


± 


0. 


11 


17.08 


± 


0. 


12 


16.99 


± 





14 


18 


17 


± 





06 


17 


07 


± 





06 


17. 


77 


± 





.06 


N5471- 


■5 


17.83 


± 





11 


17.28 


± 


0. 


12 


17.27 


± 





14 


18 


.58 


± 





06 


17 


59 


± 





06 


18. 


16 


± 





,06 


N5471- 


■6 


17.34 


± 





11 


17.21 


± 


0. 


12 


17.00 


± 


0. 


14 


17 


.62 


± 





06 


16 


41 


± 





06 


17. 


12 


± 





,06 


N5471- 


■7 


16.98 


± 


0. 


11 


16.59 


± 


0. 


12 


16.47 


± 





14 


17 


56 


± 





06 


16 


19 


± 





06 


17. 


10 


± 





.06 


N5471- 


■8 


17.89 


± 





11 


17.41 


± 


0. 


12 


17.39 


± 





14 


18 


.70 


± 





06 


17 


68 


± 





06 


18. 


22 


± 





,06 


N5471- 


■9 


19.15 


± 





13 


18.63 


± 


0. 


13 


18.63 


± 


0. 


15 


20 


.23 


± 





06 


18 


71 


± 





06 


19 


54 


± 





,06 


N5471- 


■10 


18.14 


± 


0. 


11 


17.80 


± 


0. 


12 


17.57 


± 





15 


18. 


.63 


± 





06 


17 


45 


± 





06 


18. 


22 


± 





.06 


N5471- 


■11 


18.55 


± 





12 


18.00 


± 


0. 


12 


18.11 


± 





15 


19 


.27 


± 





06 


18 


40 


± 





06 


18. 


81 


± 





,06 


NGC 5471 


14.14 


± 


0. 


10 


13.84 


± 


0. 


12 


13.74 


± 





14 


14 


.67 


± 





06 


13 


62 


± 





06 


14. 


27 


± 





,06 



Note. — Integrated magnitudes and Hct flux for the eleven regions of NGC 5471. Magnitudes from 
HST are given in the VEGAMAG system. Results for m 67 5^/ column were measured from the emission 
line free F675W image. The last row belongs to the integrated magnitud of the whole region. 



TABLE 5 
NEBULAR PHOTOMETRY 



Aperture 


Re// 


F(Ha) 


L(Ha) 


Q(H°) 


MHn 




pc 


10 -13 erg s — 1 cm~ 2 


10 39 erg s' 1 


10 51 s- 1 


10 4 M 



NGC5471-1 


61.1 


6.63 


± 





31 


1.1 


1.1 


13.0 


NGC5471-2 


60.9 


3.07 


± 





,14 


1.9 


1.9 


6.1 


NGC5471-3 


47.0 


0.93 


± 





.04 


0.6 


0.6 


1.9 


NGC5471-4 


51.2 


1.59 


± 





.07 


1.0 


1.0 


3.2 


NGC5471-5 


46.0 


0.79 


± 





,04 


0.5 


0.5 


1.6 


NGC5471-6 


58.9 


2.42 


± 





11 


1.5 


1.5 


4.8 


NGC5471-7 


63.1 


2.34 


± 





11 


1.5 


1.4 


4.7 


NGC5471-8 


43.5 


0.72 


± 





,03 


0.5 


0.4 


1.4 


NGC5471-9 


34.0 


0.38 


± 





02 


0.2 


0.2 


0.8 


NGC5471-10 


46.2 


1.22 


± 





,06 


0.8 


0.8 


2.1 


NGC5471-11 


53.6 


0.31 


± 





01 


0.2 


0.2 


0.6 


NGC 5471 


525 


36.5 


± 


1 


70 


23 


22 


73 



Note. — Integrated Ha flux and derived parameters for the eleven regions of 
NGC 5471. R e // is the effective radius of the polygon. The last row belongs to the 
integrated values of the whole region. 



We have used FUV and NUV filter passbands to cal- 
culate the AB magnitudes of Vega by means of equation 
2. As the zero points for GALEX arc ZP = 0.000 in both 
bands, we have 

m FUV = AB FUV - 2.073 (3) 

m NUV = AB NUV - 1.647 (4) 

Finally, the GALEX VEGAMAG magnitudes of 
NGC 5471 used below are rapuv = 12.12 and m^uv 
= 12.59. 

Although a complete analysis of the stellar populations 
is better characterized by the CMD study, integrated 
photometry provides additional information from wave- 
lengths where CMD studies are not feasible. However, 
the usual methodology to infer properties from clusters 
based on their integrated light is not trivial for low mass 
clusters and we need to be sure that the precedure is 
valid before the use of any analysis tool. The first triv- 
ial test is that cluster luminosity fulfills the bolometric 
assumptions implicit in the synthetic photometry (clus- 
ter luminosities larger than the Lowest Luminosity Limit 
(LLL), Ccrvino & Luridiana 2004). This basically es- 



tablishes when the cluster luminosity in a given wave- 
length/filter may be produced by just an individual star 8 . 
The mean luminosity of a cluster obtained by population 
synthesis models can be safely used when the cluster is 
brighter than 10 times the LLL, or 2.5 in magnitude units 
(Cervino & Luridiana 2004). 



The results of the LLL analysis are shown in figure 8, 
where the ages are just first order estimations when LLL 
and additional limitations in the analysis are not taken 
into account. The figure shows that, with the excep- 
tion of a few clusters, most of the clusters do not fulfill 
the 10 x LLL safety limit. Thus, formally, the analysis 
based on the mean values obtained from synthesis mod- 
els should be considered as preliminary and used with a 
cautionary note, as some could be biased by contribu- 
tions from individual stars. 

With this limitation in mind and taking into account 
that the results could be reinterpreted we have ana- 

8 The LLL just establishes this possibility, but it does not provide 
more information about whether it is really the case. 
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6.0 " 6.5 7.0 7.5 8.0 



log(Age) 

Fig. 8. — Lowest-luminosity limit for the optical and infrared 
niters used in the cluster photometry. Each cluster is located at 
the age resulting from the CHORIZOS analysis. The continuous 
line indicates de LLL and the dotted line the "safety" limit of 
LLL+2.5. 

lyzed the integrated photometric measurements with the 
help of CHORIZOS 9 , a code written in IDL by Mai'z 
Apcllaniz (2004). This code compares photometric data 
with model spectral energy distributions (SEDs), calcu- 
lating the likelihood for the full parameter range speci- 
fied, allowing for the identification of multiple solutions. 
CHORIZOS can handle SED families for stellar atmo- 
spheres and the mean SED obtained by synthesis mod- 
els 10 (Starburst99; Leitherer et al. 1999) and uses % 2 
minimization to find all models compatible with the ob- 
served data in a (2 intrinsic +2 extrinsic)-dimensional 
parameter space. It includes a wide list of filters for dif- 
ferent systems and allows for published magnitudes to 
be used directly. For a set of objects with M + 1 mea- 
sured magnitudes (with independent uncertainties), de- 
rived M independent colors, and a set of N parameters 
with M > N, no strict solutions, but only approximate 
ones, can be found since the problem has more equations 
or sets of model parameters that produce a given color 
than model parameters. Hence, the code finds the N- 
volume of approximate solutions by using x 2 minimiza- 
tion by the comparison with synthetic photometry. In 
the case of stellar clusters the synthetic photometry is 
obtained from the average SED of population synthesis 
models. In these cases, and in addition to the success of 
the 10 x LLL test, the x 2 minimization procedure is valid 
as far as (a) the synthetic photometry used has an in- 
trinsic variance due to size of sample effects much lower 
than the uncertainty in the data (Cervino & Luridiana 
2006), (b) size of sample effects do not produce biased 
synthetic photometry (Cervino & Valls-Gabaud 2003), 

The latest version can be obtained from 
http : / /www. iaa.es/ ^jmaiz 

10 Synthesis models provide an average SED with a relative vari- 
ance equal to zero only in the case of clusters with an infinite num- 
ber of stars Cervino & Luridiana (2006). 



and (c) the average values obtained from the models and 
used in the fit are representative of the intrinsic under- 
lying distribution, that is, the underlying distribution 
of integrated luminosity have a gaussian-likc shape and 
the intrinsic variance can be translated to gaussian-likc 
confidence intervals and x 2 minimization procedures are 
valid (Cervino & Luridiana 2006) 11 ; this last require- 
ment implies the analysis of high-order moments of the 
distribution and a first order estimation of the mass and 
age of the cluster in absence of external constrains and 
once the 10 x LLL requirement is fulfilled by the observa- 
tional data. Although, from the LLL test, we know that 
the results from CHORIZOS should be interpreted with 
caution for the individual clusters analyzed (but valid 
for the cluster as a whole), we have used it to obtain 
a first order estimation of the effects of extinction and 
reddening which would modify the LLL test. 

The first step of the program is to read the unreddened 
SED models, extinguish them according to a given ex- 
tinction law, and obtain the synthetic photometry. To 
run the program, CHORIZOS requires a file with the 
photometry data and a series of input parameters, such 
as the model family, understood as the collection of SED 
from stellar atmosphere models or synthetic stellar pop- 
ulations, and the range of the parameters to be used. 
Then, CHORIZOS reads the model photometry from the 
tables calculated previously, interpolates, calculates the 
likelihood for each grid point, and writes the results for 
each of the four parameters. 

In the case of our study, we have used the STAR- 
BURST99 models (Leitherer et al. 1999) with the im- 
plementation by Smith et al. (2002) which includes a 
grid of blanketed stellar atmospheres at different metal- 
licities. We fixed two parameters: the known metallicity 
of NGC 5471 and the type of dust, leaving unconstrained 
the amount of extinction and the age. 

Regarding extinction, CHORIZOS works with 
i?(4405 — 5495) and -R5495, which are the monochromatic 
equivalent to E(B — V) and Ry respectively, the values 
4405 and 5495 referring to the central wavelength (in 
A) of the B and V filters. We used the standard type 
of dust characterized by a value .R5495 = 3.1. The 
metallicity of the region was assumed to be traced 
by the oxygen abundance. Recent derivations from 
detailed spectrophotometric studies give a mean value 
of 12+log(0/H) of 8.06 ± 0.03 for components A, C and 
D; a slightly lower abundance of 7.94 ± 0.03 is found for 
component B (Kennicutt, Bresolin & Garnett 2003). 
These values correspond to about 0.25 the solar oxygen 
abundance, if the value by Asplund et al. (2005) for the 
solar oxygen abundance is assumed. Therefore we used 
the value Z = 0.2Z Q in CHORIZOS. 

4. DISCUSSION 

In this section we discuss the results of the photometry 
in terms of the resolved stellar photometry in the CMD 
and in terms of the integrated properties, both of the in- 
dividual clusters defined in section 3.2 and of the whole 
star forming complex. We attempt to bring all the re- 
sults into a coherent picture of the past and present star 
formation processes in NGC 5471. 

This requirement is actually linked with the lOxLLL and the 
unbiased photometry requirement. 
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4.1. Cluster Analysis 

Analysing the integrated properties of a complex re- 
gion like NGC 5471 is an exercise that has to be un- 
dertaken with care, because the integrated photometry 
includes a wide range of regions with different cluster 
masses, ages and extinction values. Thus the result of 
synthesizing the integrated ultraviolct-to-NIR photome- 
try with CHORIZOS, which finds the most likely single 
stellar population (SSP) under unconstrained inputs for 
mass, age and extinction, is only of general interest as 
a guiding signpost. If we use only the optical-to-NIR 
photometry, CHORIZOS finds a most likely SSP with 
relatively high extinction and a young age 3 ± 2 Myr. 
If we include the FUV and NUV GALEX photometry 
then CHORIZOS finds two loci of maximum likelihood 
for a single SSP, one that corresponds to an age of 8 ± 2 
Myr with moderate extinction and a second maximum at 
an age 60 ± 15 Myr and very low extinction. This result 
suggests that a significant fraction of the GALEX flux in- 
tegrated within the 15 arcsec radius aperture from NGC 
5471 is produced by a post-nebular stellar population of 
B stars (see also Waller et al. 1997). 

Three previous studies have modelled the age of NGC 
5471 using ultraviolet-to-optical spectra. Rosa & Ben- 
venuti (1994) use an HST/FOS spectrum through an 
aperture of 4.3 arcsec around knot A. Mas-Hesse & 
Kunth (1999) use an IUE spectrum through an aper- 
ture of 10x20 arcsec and the FOS optical spectrum of 
Rosa & Benvenuti (1994). Pellerin (2006) uses a FUSE 
spectrum through an aperture of 4x20 arcsec. All three 
works find that the UV spectrum of NGC 5471 can be 
explained with a young burst of less than about 5 Myr 
and a relatively low extinction. These results are consis- 
tent with the detection of Wolf-Rayct (WR) stars in NGC 
5471 (Mas-Hesse & Kunth 1999). However these previous 
works may not be sensitive to an older population since 
they are based on the analysis of UV continuum and ab- 
sorption lines (Si iv and C iv) and optical emission lines 
as primary age indicators, and are thus biased towards 
younger knots with brightest ultaviolet fluxes; while the 
integrated GALEX photometry through a large aperture 
is more sensitive to the extended, fainter, more evolved 
stellar population 12 . Indeed, Meurer et al. (1995), Maoz 
et al. (1996), and Chandar et al. (2005) have found that 
the relative contribution from the cluster stellar popula- 
tion to the total UV flux in starbursts is only a fraction in 
the range -20% to —40%, the majority of the UV flux 
coming from a diffuse fainter stellar population that is 
most prominent when the total flux is integrated, just as 
in the case of NGC 5471. In addition, the analysis per- 
formed with CHORIZOS focus the results on the contin- 
uum shape. This discrepancy between ages inferred from 
emission lines and from the continuum points towards a 
complex star formation history in NGC 5471 as we show 
below. 

The integrated photometry analysis of the clusters 
identified in section 3.2 may be affected by size of sam- 
ple effects due to the failure of the lOxLLL test in some 
of the clusters. However, we present here the analysis 
to show, as a case of study, the possible interpretations. 
The values quoted must not be taken, in any case, as 

12 Recall that the la detection for the GALEX data is 
ABpuv=27.6 mag arcsec -2 
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Fig. 9. — Ho and J-V versus J-H plot for the clusters identified 
in section 3.2 (V stands for the 547M filter). A clear correlation of 
decreasing Ha with redder colors points to an aging effect within 
the ongoing starburst phase. 

final values for cluster ages. 



Figure 9 shows a color-color-Ha plot of the individual 
clusters. The hypernova knot is not shown, and knots 3 
and 4 have been combined into a sing lc knot 34 13 . There 
is a clear correlation in the sense that redder knots have 
less Ha flux. The flux of Ha can decrease due to the 
aging of the burst or to different clusters having differ- 
ent masses, less massive bursts producing fewer ionizing 
photons 14 . This latter could be the case, including a re- 
lationship with color, when strong size of sample effects 
are present and just a few stars are responsible for most 
of both the Ha and continuum emission (over-luminous 
case) or by the absence of the most massive stars which 
mimics variations in the upper mass limit of the Initial 
Mass Function (under-luminous case; e.g. Villaverde et 
al. 2010). 

Knot 1 is the brightest in Ha and the age estimations 
obtained by CHORIZOS give the youngest age, 3 ± 2 
Myr, with the lowest extinction. For knots 2, 3, 4, 7, 
8, 9, 10 CHORIZOS estimates an intermediate age, 5 ± 
2 Myr, with an intermediate extinction, except may be 
knot 9 which has some indications for a somewhat higher 
extinction. Knots 5 and 11 have older age estimations, 
9 ± 3 Myr. However, except for Knot 1, 2 and 34, the 
age estimations can also be explained by size-of-sample 
effects with under-luminous case. In these cases, age 
estimations, although apparently consistent, should not 
be considered as final, pending a more complete analysis 
with additional data sets. 

13 The nature of knots 3+4 is not clear, it may be that both 
are part of the same cluster that has evolved sufficiently for the 
stellar winds to evacuate its center from ionized gas. Further higher 
spatial resolution IR and UV imaging is needed to establish more 
details. 

14 Extinction may also play a small role but values significantly 
larger than reported in the literature, which are lower than E(B- 
V)=0.14, would be required. 
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In summary, from the cluster analysis of individual 
bright knots and from the integrated photometry of the 
whole complex there are clear indications that there is a 
complex star formation in NGC 5471 event, with indica- 
tions of an older population of B stars and a luminous 
young ionizing cluster. This will be well established with 
the analysis of the CMD in the next section. 

4.1.1. Mass of the ionized gas 

The H + mass and the mass of the ionizing cluster may 
be derived from the Ha flux. The published values of in- 
ternal extinction for NGC 5471 differ from different au- 
thors, but they are lower than about E(B-V)=0.14 on av- 
erage (Skillman 1985; Torres-Peimbert et al. 1989; Ken- 
nicutt & Garnett 1996; Luridiana et al. 2002; Esteban et 
al. 2002). In NGC 604 Mafz Apellaniz et al. (2004) have 
estimated that only about 10% of the Ha flux is miss- 
ing in the higher opacity zones. Therefore, we have not 
attempted to correct the measured Ha flux to estimate 
other magnitudes from this value. The electron density 
asummcd is n e = 100 cm" 3 (Luridiana et al. 2002). 

The Ha flux in a 15" radius circular aperture (uncor- 
rected for internal extinction) yields a lower limit to the 
luminosity Ln a ~ 2.3xl0 40 ergs" 1 . The total number of 
H-ionizing photons, Q(H°), estimated from Ostcrbrock 
(1989) is: 

Q(H°) = ^jx P^- ~ 2.2 x lO 52 *- 1 

ri a 

where Q-bI^hL ~ 2.99, assuming case B recombination, 
olb ~ 2.59 x 10~ 13 cm -3 s _1 , and an effective recombi- 
nation coefficient for Ha of 8.63 x 10 -14 cm -3 s _1 for 
10 4 K and n e = 100 cm" 3 (Osterbrock 1989). Thus the 
total (lower limit) ionized mass of the gas cloud is 

M H + = Q(H )-^- ~ 7.3 x 10 5 M S 
n e aB 

From evolutionary models of single stellar population 
ionizing clusters and radiation-bounded H II regions, a re- 
lation between Q(H°) per solar mass and the H/3 equiv- 
alent width can be found (Diaz 1998), which allows to 
estimate the stellar cluster mass by taking into account 
the cluster evolution: 

log[Q{H°) /M Q ] = 0.86 log[EW(HP)] + 44.48 (5) 

Hence, a lower limit for the mass of the cluster can 
be estimated by means of the H/3 equivalent width and 
the Ha luminosity of the region. Using the range of val- 
ues EW(H^)= 134-296 A published by different authors 
(Luridiana et al. 2002; Esteban et al. 2002; Rosa & Ben- 
venuti 1994; Rosa 1981), and the value of Q(H°) cal- 
culated above, equation 5 yields a mass of the ionizing 
stellar cluster in the range 0.6 — 1.1 x 10 6 M©. Using 
the calibration of Kennicutt (1998) and the luminosity 
that we measure, Ln a ~ 2.3 x 10 40 erg s _1 , we obtain 
a star formation rate of 0.18 M© yr" 1 , or 2.1 x 10~ 7 
Mq yr -1 pc -2 (taking the area of the region as the cir- 
cular aperture defined in Figure 4), that translates into 
a cluster mass of 10 6 M Q for a burst duration of 10 Myr, 
which is the mean duration of the current burst implied 
by the analysis in the previous section. Nevertheless, it 



is known that this calibration was computed for a partic- 
ular star formation episode at nearly constant rate dur- 
ing more than 100 Myr, so it might be not suitable for 
bursty systems. We have applied the new calibration by 
Otf-Floranes & Mas- Hesse (2010) 1 , which makes a dis- 
tinction between instantaneous (IB) and extended burst 
(EB). We used as input our measured Ha luminosity, no 
color excess, and a Salpeter Initial Mass Function in the 
range 1-100 M Q . For the case of instantaneous burst, 
Oti-Floranes & Mas-Hesse (2010) use the concept of star 
formation strength (SFS), defined as the total mass of 
gas transformed into stars during the burst, giving three 
different ages for the IB case. We obtained a SFS of 1.5 
x 10 6 M for a burst of 4 Myr, 2.3 x 10 6 M Q for 5 Myr, 
and 1.0 x 10 7 M Q for the 6 Myr case. For a EB, the 
results are given in terms of the stardard star formation 
rate (SFR), with three different tabulated ages, namely: 
10, 30 and 250 Myr. However, in this age range the re- 
sults are independent of the age for the Ha estimator. If 
the region would be experiencing a continuous constant 
star formation, the SFR given by this calibration would 
be 0.1 M© yr -1 , less than a factor of 2 lower than the 
calibration by Kennicutt (1998) 16 . 

Thus we find that the masses of ionized gas and of 
the stellar cluster in the current star formation event 
arc both of the same order of magnitud, close to 10 6 
M©. This similarity may be understood in terms of a 
very high efficiency in star formation. Given that this 
is not the first event of star formation in NGC 5471 in 
the last 100 Myr, this high star formation efficiency takes 
an important meaning that implies that there has to be 
a large reservoir of gas to maintain this important star 
formation activity. 



Although optical images show NGC 5471 apparently 
isolated in the outskirts of M101, the Hi distribution of 
M101 shows that NGC 5471 is, in fact, one of the bright- 
est regions in Hi, located in a spiral arm (cf. figure 1). 
The thermal radio continuum emission from NGC 5471 
has the highest observed central emission measure of all 
the M101 complexes. Closely coincident with this Hn 
radio continuum peak, there is an Hi concentration of 
nearly 10 8 M (Viallefond et al. 1981). The general dis- 
tribution of the H I gas in M101 follows that of the diffuse 
FUV emission, with enhancements in the 21 cm line emis- 
sion seen near FUV peaks. Smith et al. (2000) include 
NGC 5471 as one of their 35 photodissociation regions 
(PDRs) candidates. Their Ultraviolet Imaging Telescope 
(UIT) FUV flux measurement is in good agreement with 
our GALEX measurement of 9.8 x 10 -14 erg cm" 2 s" 1 
A -1 , and NGC 5471 is clearly the brightest FUV peak 
in the whole M101 galaxy. This is reflected in a de- 
creased molecular content surrounding NGC 5471 which 
has been dissociated by the high star formation efficiency, 
the highest in M101 as calculated by Giannakopoulou- 
Creighton et al. (1999). 

4.2. CMD 

15 Their web tool can be found at 
http : //www. laef f . cab . inta-csic . es/research/sf r/ 

16 This calibration integrates a Salpeter IMF in the mass range 
0.1-100 Mq, which predicts a larger mass. 
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Fig. 10. — Final CMD with isochroncs divided into five age bins, drawn by their concave hull, and shown with different colours. Error 
bars correspond to an average error in each magnitude bin. The right hand side of the figure shows the luminosity function. The line across 
the bottom shows the completeness limit at the 50% level in the core (dashed line) and in the rest of the image (solid line). 



The photometry of the individual stars presented in 
the color-magnitude diagram, vci^ym— ni675W (V-R) vs. 
111547M (V), of figure 10 is a powerful tool to further un- 
derstand the star formation history in NGC 5471. The 
isochrones represented in the figure have been divided 
into five age bins, (4-10, 10-25, 25-50, 50-75, 75-100) 
Myr, and drawn by their concave hull, shown with differ- 
ent colours. The right hand side of the figure shows the 
luminosity function, which peaks at M547A/ = 24, and 
becomes incomplete for fainter stars, as indicated by the 
50% completeness limit shown. We recall that this CMD 
is the result of a careful process of calibration, measure- 
ment and cleaning, as detailed in section 3.1 above. We 
have not included in the CMD those point-like sources 
with a partially resolved PSF, that are likely to be com- 
pact clusters. At the distance of M101 some of the point 
sources are unresolved multiple systems, some of which 
may have been included in the CMD, but these will not 
really affect any of the arguments that we present here. 

A glance at the CMD tell us at once that star forma- 
tion has been proceeding more or less continuously for 
the last 100 Myr. The diagram is populated by stars of 
^4 — 50 Mq. We can readily identify clear structures 
in the diagram. A well-defined main-sequence at V-R 
~ -0.1, which begins with the brightest star at 19 mag 
(that corresponds to an age of ~ 4 Myr) and goes all the 



way down to 25 mag. Red stars span a range of 3 — 4 
magnitudes in brightness, indicating that star formation 
has occurred during the last ^15 — 70 Myr. There are 
some traces of the blue-plume at My = —9, indicating 
ongoing star formation. A strip of stars following a patch 
from [0.6,25] to [1.1,22.5] form the red giant branch. 

4.2.1. Mass in stars and star formation history 

In order to quantify the relative importance of the 
present burst to the star formation history during the 
previous 100 Myr in NGC 5471, we make a rough esti- 
mate of the stellar mass from the CMD. We follow three 
different approaches. One method consists in assigning 
a mass to each of the stars in the CMD, corresponding 
to the most probable isochrone. This is done considering 
a 2D Gaussian probability density distribution centered 
on the distance in the color-magnitude space of a given 
observed star to the isochrones with dispersions being 
the estimated errors for the photometry. Once a mass 
has been assigned to a star, we count the stars in the 
mass range 6-55 M and integrate a Salpeter IMF in 
the mass range 1-100 M Q . The result of this calculation 
gives ~ 2 x 10 5 M©. The second method is as follows: 
first we assign a mean stellar mass of 6 M Q to the stars 
in the most populated bin of the luminosity function, 
111547M = 23.5. This mass is then scaled up by the ratio 
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Fig. 11.— STARFISH recovered Star Formation History of NGC 
5471. 

of the total 111547^/ of NGC 5471 to the magnitude of the 
bin ni547M = 23.5. The resulting mass is used to nor- 
malize the IMF at the bin corresponding to 6 M Q . The 
result of integrating a Salpeter IMF in the range 1-100 
M© thus normalized gives a value of ~ 5 x 10 5 M©. 



A third independent method to estimate the mass is 
based on the fit of the CMD by means of Hess diagrams. 
We have used the code STARFISH 17 (Harris & Zaritsky 
2001, 2002), with the lastest patchsets, to find a best fit 
to the CMD by a linear combination of Hess diagrams at 
different ages. We have used our own set of isochrones de- 
scribed in section 3.1.3. The best-fit SFH of NGC 5471 is 
shown in Figure 11. Figure 12 shows the synthetic CMD 
of NGC 5471 constructed using the SFH amplitudes out- 
put of the best-fit SFH; the observed CMD is also shown 
for reference. The synthetic main sequence seems tighter 
and bluer compared to the observed one, probably due 
to differential reddening and underestimated errors. The 
SFH derived seems to be consistent with a continuous 
one, but the errors only allow to take into consideration 
two episodes: a dominant young burst < 20 Myr, pre- 
ceded by a longer event during the period 20-100 Myr. 
According to the SFR, the integrated mass in stars of 
the youngest event is comparable to the oldest one, with 
~ 10 6 M . 

Although these methods give only an order of magni- 
tude estimate, we can conclude that the integrated mass 
of stars formed in the past 20-100 Myr is of the order of 
magnitude of the current star formation event, that we 
have calculated above as ~ 10 6 M©. If the star formation 
proceeds in an isolated gas cloud, we would expect that 
the rate at which stars form would decrease with time, as 
less gas is left available to form new stars. This is what 
is observed in many giant H n regions with two or more 
generations of star formation, where the latest event is 
less massive than the previous. In the case of NGC 5471 
the results from the integrated photometry and from the 
CMD indicate that the current star forming event is as 
massive as those during the previous 20-100 Myr. We can 

17 Available at http://www.noao.edu/staff/jharris/SFH/ 
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Fig. 12.— NGC 5471 synthetic CMD based on the best-fit SFH 
solution of STARFISH (left) and the observed one (right). 

understand this in the context of the large H I spiral arm 
in which the region is immersed (cf. section 4.1.1). In 
this respect, Waller et al. (1997) conclude from a morpho- 
logical study of M101 that NGC 5471 may be the result 
of tidal interactions of M101 with the nearby galaxies 
NGC 5477 and NGC 5474 in the time interval 100-1000 
Myr ago. 



4.2.2. Spatial distribution of the star formation history 

Given the long duration of the star formation in NGC 
5471, ~ 100 Myr, and its large size ~ 1 kpc, it is natu- 
ral to question the spatial location of the star formation 
through time along the extent of the region. In order 
to visualize this, we have selected four different areas in 
the CMD and proceeded to locate spatially the stars in 
these four boxes. Figure 13 displays the resulting distri- 
butions. The youngest and more massive stars (top-left 
panel) are clearly concentrated in the main star forming 
emission line knots, implying that massive stars are now 
forming mainly in clusters. Only two stars are not di- 
rectly associated with these knots, one of which belongs 
to an ishocronc of 20 Myr. The bottom-right panel dis- 
plays the distribution of intermediate mass stars (~ 5—10 
M©) older than about 50 Myr; these stars are distributed 
mainly towards the halo, with only a few in the core of 
the region. The top-right and bottom-left panels show 
the distribution of intermediate mass stars (~ 10 — 15 
Mq), both those which are in the main sequence and 
those which have already evolved out of it. The distri- 
butions of stars in these two CMD boxes are clearly con- 
centrated fairly uniformly throughout the core of NGC 
5471; some of the stars in these two boxes are found in 
the halo, with a sparse distribution somewhat similar to 
the older lower mass stars in the bottom-right panel. Be- 
cause the stars in these two intermediate boxes arc of the 
same mass but they are found in a range of evolutionary 
stages, this implies that star formation has proceeded 
more or less uniformly in the core for the last ^20 — 50 
Myr. 

From this analysis of the spatial distribution of the 
CMD, a clear picture emerges in which the star forma- 
tion in the complex NGC 5471 has proceeded in a general 
spatio-temporal sequence from the periphery to the core. 
During the first epoch, ~ 50 — 100 Myr ago, the star for- 
mation occurred along the periphery and the core, then 
a second major event ~ 10 — 20 Myr ago was more or 
less uniformly distributed mainly in the core, and finally 
the current ionizing star forming event is mostly concen- 
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Fig. 13. — Spatial distribution of the star formation with time. The CMD contains only stars with magnitudes below the 50% completeness 
defined in the core region. The spatial distribution of the stars within four boxes selected in the CMD are shown in the right hand side 
images. The images are also drawn with a circle of 15 arcsec radius, indicating the area of the integrated photometry, an ellipse delineating 
the core region where the present star forming event occurs, and the polygons selected for the cluster analysis. The square points represent 
the stars in the box selected in the CMD, the light color points indicate the stars located in the core of the region (inner ellipse), while the 
black points stand for the stars belonging to the halo. Notice how the younger stars are concentrated in the inner clusters (top-left panel) 
and intermediate and older stars populate uniformily the region. 



tratcd in the singular bright clusters located inside the 
core, where the gas reservoir is presumably higher 18 . At 
the same time, we have seen from both the integrated 
photometry analysis and the CMD, that as the events 
proceeded concetrating from the periphery to the core 
clusters, the amount of mass of gas formed into stars 
seems to be of the same order of magnitude, but more 
concentrated in time, with higher star formation rate. 

We note that the current event of ionizing clusters is 
apparently contained within a large bubble which defines 
the core of NGC 5471. This bubble, of projected size 
400 x 550 pc, is clearly defined geometrically (cf. figure 
2) and kincmatically (Muhoz-Tuhon et al. 1995), and is 
likely to have been produced by the stars that formed 
~ 20 Myr ago. 19 

There are other examples that seem to follow the gen- 
eral process of star formation that we have outlined here 
for NGC 5471. Ubeda et al. (2007) find an inwards se- 
quence of star formation and a bubble for cluster I-A in 
NGC 4214. Walborn & Blades (1997) find a triggered 
spatio-temporal distribution of stars in 30 Doradus, but 
this sequence occurs at the scale of the individual clus- 
ters in NGC 5471. Indeed there are indications that at 

18 Notice that the age dating results from the individual clus- 
ter photometry is compatible with this age-spatial distribution of 
the younger stars in the CMD, and this reinforces the age results 
obtained for the clusters. 

19 The results from a kinematical analysis of the expansion of 
this bubble arc in progress and will be published elsewhere 



this smaller scale, ^50pc, the evolution of the different 
individual clusters indicated by the age dating is also re- 
flected in their morphology (e.g. as mentioned in the case 
of 3+4) , but a detailed study at these scales requires the 
high resolution UV-optical-IR imaging that can only be 
provided by the WF3 camera onboard HST. 

5. CONCLUSIONS 

We have performed a photometric study of NGC 5471, 
a GEHR in the outskirts of the spiral galaxy M101, to de- 
rive its star formation history by means of resolved stellar 
photometry and integrated cluster photometry. Further- 
more, integrated photometry integrated photometry of 
the whole region, using data from GALEX (ultraviolet), 
HST/WFPC2 (optical), and TNG (near infrared), yields 
two possible solutions for a single stellar population: one 
that corresponds to an age of ~8 Myr with moderate ex- 
tinction and one with an age around 60 Myr and very 
low extinction. This degeneration is typical of this kind 
of approach when spatially integrated data are used and 
reflects the complexity of GEHR which include clusters 
with different ages and extinction values. 

From the photometric analysis of the eleven clusters 
defined on the IR H image a correlation emerges in the 
sense that redder knots have less Ha flux, showing a clear 
aging trend. The ages range from 3 Myr for the youngest 
cluster, up to 10 Myr for the oldest one. However, the 
Lowest Luminosity Limit test prevents us from taking 
these estimates at face value, since only a few individual 
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clusters satisfy the required limit. 

The complex history of star formation of NGC 5471 
revealed by the cluster analysis is confirmed by the re- 
solved stellar analysis. From the CMD it is clear that 
star formation has been proceeding continuously with 
more or less intensity for the last 100 Myr. The well- 
defined main-sequence gives a youngest age of about 4 
Myr, while red stars, which span a range of 3-4 magni- 
tudes, indicate that star formation has occurred during 
the last 15-70 Myr. 

We have found that the masses of the ionized gas and 
of the stellar cluster in the current star formation event 
are both of the same order of magnitude, around 10 6 
Mq, which may be understood in terms of a very high 
efficiency. In fact, the star formation strength of NGC 
5471 is more than a factor of 2 that of NGC 604 in M33 
(Relano & Kcnnicutt 2009). This implies that there has 
to be a large reservoir of gas to maintain this star for- 
mation rate, something that is confirmed by the huge Hi 
concentration in which the region is located. 

From the spatial distribution of the stars, we can con- 
clude that the star formation in NGC 5471 has proceeded 
in a general spatio-temporal sequence which was more 
extended in the past and now is concentrated in clus- 
ters in the core. The current event of ionizing clusters is 
contained within a large bubble, which is likely to have 
been produced by the stars that formed <~ 20 Myr ago. 
This leads to the general question of wether a low density 
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"seed" generation of evolved intermediate mass stars is in 
fact a prerequisite to trigger a starburst event (Brandner 
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lution, such as that can be provided by the WF3 camera 
on board HST, are obviously needed in order to have a 
conclusive answer. 
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